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ABSTRACT 



Aims. The population of solitary compact objects in the Galaxy is very diffcult to investigate. In this paper we analyze the possibility 
of using microlensing searches to detect and to analyze the properties of the solitary black holes and neutron stars. 
Methods. Evolution of single and binary stars is considered using the StarTrack population synthesis code. We investigate the prop- 
erties of the Galactic population of compact objects numerically. 

Results. We find that the compact object lensing events are concentrated in a region with the radius of « 5 degrees around the 
Galactic center. The distribution of masses of the lenses for the models we consider differs but only slightly from the underlying mass 
distribution. The expected detection rates are of the order of a few per year. 
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1. Introduction 

The stellar evolution models predict the Galactic population of 
compact objects - neutron stars and black holes - numbers some- 
where between 10^ and 10'' objects. Most of them are solitary 
and only a small fraction resides in binaries. A small fraction 
of solitary neutron stars is visible as radio pulsars. Some black 
holes in binaries are actively accreting and are observable mainly 
as X-ray transients, however most of the black hole population 
has not been yet seen. The properties of black holes population 
in our galaxy depend on the history of star formation rate, evolu- 
tion of metallicity and on the initial mass function and on details 
of compact object formation in supernovae explosions. 

The properties of the population of solitary black holes 
can only be investigated indirectly, through observations of 
their interaction with the interstellar matter or light emitted by 
stars. Solitary black holes should be accreting gas from the 
ISM. Therefore some of them should be observable i n X-rays 
jAgol & Kamionkowskill2002l:lBeskin & Karpovll2005h . The lu- 
minosity in X-rays shall on one hand be smaller than for solitary 
neutron stars where the surface emission plays a significant role, 
but on the other hand it may be increased due to smaller veloci- 
ties and higher mass of black holes in comparison with neutron 
stars. The searches for such objects have not been successful yet. 

Solitary black holes are also d etectable wi th current mi- 
crolensing s earches(lPaczvnskil2003h like OGLE (IUdalskil2003h 
and MOA jBond et al.l 2001b . These campaigns have already 
yielde d several candidate black holes detections. Bennett et al.l 
I2OO2I) has presented two events: MACHO-96-BLG-5 and 
MACHO-98-BLG-6, with the mass estimates 6'^1°Mq and 
6^1Mq respective! v. iMao et al.l(l2002l) showed that OGLE-1999- 
BUL-32 identified also as MACHO-96-BLG-22 is a black 
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hole candidate with the minimum mass of lO.SMp. A further 
search for X-rays from MACHO-96-BLG-5 (iMaeda et al.ll2005l) 
yielded an upper limit corresponding to the luminosity of less 
than (8 - 9) X 10'^ ° LEdd- A recent l ikelih ood analysis of 22 
microlensing events Poinde xter et al.l (l2005l) lead to confirma- 
tion of the black hole candidate MACHO-99-BLG-22, while the 
other candidates are less probable. 

The sample of the microlensing compact objects and black 
hole candidates will increase with time. It is therefore interesting 
to see what constraints can be imposed on the models of black 
hole formation and evolution by these observations. In this pa- 
per we present a simulation of the stellar evolution leading to 
production of solitary black holes in our galaxy. We analyze two 
scenarios: the single stellar evolution and the fromation of soli- 
tary black holes through disruption of binaries. We then analyze 
the motion of black holes in the galactic potential and search for 
possible microlensing observable from the Earth. 

In section 2 we present the basic ingredients of the calcula- 
tion: the stellar evolution model, the galactic potential and mass 
distribution used, and the lensing search algorithm. Section 3 
contains the results and section 4 the discussion. 



2. Description of the model 

2.1. Compact object formation 

The single star ev olution is modeled using the formulae of 
iHurlev et al.l (l2000l) and we use the S tarTrack population syn- 
thesis code dBelczynski et al.l[2002ah for modelling binary evo- 
lution. The StarTrack population synthesis code was initially 
developed for the study of double compact object mergers in 
the context of GRB progenitors (Belczynski et al. 2002aj^and 
gravitational-wave inspiral sources dBelczynski et alj2002bl) . In 
recent years StarTrack has undergone major updates and re- 
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Table 1. Models of black hole formation considered in this paper 

Model Description 

A Standard 

W05 stellar winds decreased by 0.5 

KO Black holes receive no kicks 

Kl Black holes receive full kicks 

C Hertzsprung gap stars can be donors in CE phase 

S Only single stars 

M Black holes from mergers 



visions in the physical treatment of various binary evolution 
phases. The new version has akeady been test ed against obser- 
vation s and detailed evolutionary calculations (tBelczynski et al.l 
|2006|) . and has been used in various applications. The most 
important updates for compact object formation and evolution 
include: a full numerical approach to binary evolution due to 
tidal interactions and coupling calibrated using high mass X- 
ray binaries and open cluster observations, a detailed treatment 
of mass transfer episodes fully calibrated against detailed cal- 
culations with a stellar evolution code, updated stellar winds 
for massive stars, and the latest dete rmination of natal kick 
velocity distribution for neutron stars dHobbs et al.ll2005T) . For 
disrupted (by supernova explosion) bina ries we follow trajec - 
tories of components as described in Bel czynski et al.l (|2006|) . 
In the helium star evolution, which is of a crucial importance 
for the formation of new classes of double compact objects, 
e.g. Jlvanova et al.l (l2003h . we have applied a conservative treat- 
ment matching closely the results of detailed evolutionary cal- 
culations. The NS-NS progenitors are followed and checked 
for any potential Roche lobe overflow (RLOF). While in the 
mass transfer phase, systems are examined for potential devel- 
opment of dynamical instability, in which case the systems are 
evolved through a common envelope phase. We treat common 
envelope events through the energy formalism ( Webbink 1984; 
iBelczvnski et al.l l2002bl) . where the binding energy of the en- 
velope is determined from the set of He star model s calcu lated 
with the detailed evolutionary code by llvanova et al ] (l2003h . For 
some systems we observe, as before, extra orbital decay lead- 
ing to the formation of very tight short lived double compact 
object binaries. However, since the progenitor evolution and the 
final RLOF episodes are now followed in much greater detail, 
we note significant differences from our earlier studies. For a 
detailed descript i on of the revised code we refer the reader to 
IBelczvnski eTaTI ( l2006l) . 

The simplest scenario leading to a single compact object is 
through evolution of single stars. In this case we put the newly 
formed star in the galactic disc with its velocity drawn from an 
appropriate distribution after the supernova event. 

An important additional scenario follows a disruption of a 
binary as a result of a first supernova explosion. A large frac- 
tion of binaries containing massive stars is disrupted in this way. 
This leads to a formation of a single compact object and a sin- 
gle companion which may still be massive enough to produce 
another compact object. We assume that the first supernova ex- 
plosion takes place in the galactic disk. We follow the motion of 
the single star until it produces another compact object. 

A much smaller contribution comes from the disruption of a 
binary during the second supernova explosion. In this case two 
single compact objects are formed. We follow the motion of the 
system after the first supernova until the second one. Thus the 
newly formed single compact objects start at the actual location 
of the second supernova with appropriate velocities. 



Finally, we also investigate black hole formation through 
mergers. This includes mergers of stars during their nuclear evo- 
lution, as well as mergers of a compact object with a massive 
companion as a result of binary interaction. In the first case we 
assume that a single star is formed with the mass equal to the 
sum of the masses of the two merging components and we evolve 
this star neglecting the potential changes in metallicity. In the 
latter case we assume that a single black hole is formed with 
the mass equal to the sum of the masses of the compact object 
and the core of the companion while the envelope is expelled. 
The newly formed black holes do not gain additional velocities 
during mergers. 

We ignore the contribution of double black hole systems, and 
single black holes formed in mergers of double black holes as 
they represent a much smaller population than the ones men- 
tioned above. 

We consider several models of binary evolution to asses the 
sensitivity of the results. We decrease the strength of the stel- 
lar winds which affects mainly the masses of the newly formed 
compact objects (model W05). A second parameter that may in- 
fluence the results significantly is the distribution of the kicks re- 
ceived by the newly formed compact object. Apart from the stan- 
dard model in which the value of the kick is decreased with in- 
creasing fall-back mass, we consider two extreme cases: in first 
black holes receive no kicks at all (model KO), and in the second 
one kick distribution regardless of their mass is the same as for 
neutron stars (model Kl). Finally in model C we allow for the 
common envelope evolution initiated by stars passing through 
the Hertzsprung gap. 

The list of the models considered is shown in Table [T| 

2.2. Galactic model 

We consider a model of the Galaxy consisting of three com- 
ponents: bulge, disk, and halo. The bulge and disk potential 
are described by the iMivamoto & Nagail (Il975h type potential 
(lPaczvnskillI990l:lBulik et alJI 19991) 
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^jR^ + (a+ Vz^ -I- Z72)2 



(1) 



where M is the mass of a given component, R - ^Jx^ + y^, and 
fl, B are the parameters. The halo is described by the density dis- 
tribution p = pc[l + (r/ri;)^]^^ with a cutoff at r„„ - 100 kpc 
above which the halo density is zero. The corresponding poten- 
tial for r < rcut is 
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We use the f oll owing values of the parameters 
aes & Rajagopall Il99lh describing the potential: 
ai = Okpc, bi = 0.277 kpc, 02 = 4.2 kpc, b2 = 0.198 kpc. 
Ml = 1. 12x10'° Ma, M2 = 8.78 X 10'° Mo, M,, =5.0x10'° Mq, 
and rc - 6.0 kpc. 

The distri bution of stars i n the disk is assumed to be that of 
a young disc (iPaczvnskill 199 0). The radial and vertical distribu- 
tions are independent i.e. the distributions factor out: 



P(R,z)oc R(R)dRp(z)dz, 



(3) 



where the radial distribution is p{R) oc exp(-R/Re^p), and Rg^p = 
4.5 kpc, and we introduce an upper cutoff at /?,„av = 20 kpc. 
The vertical distribution is exponential p{z) oc exp(-z/75pc) 
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Providing that A = 4(D ■ Pf - 4D^ {p^ - R^) > the Eq. Q 
has a formal solution: 



Fig. 1. Geometry of the lensing event. Star and the lens are lo- 
cated inside a sphere of radius R equal to the radius of the Earth's 
orbit. Vector P denotes the position of the star in the galaxy, 
while £) is a star-lens vector We look for a situation when star- 
lens line crosses the Earth's orbit during the evolution of the sys- 
tem. 



. This is not a self consistent approach since the matter den- 
sity corresponding to the disk potential is not the same as the 
stellar density. A self consistent approach has been investi- 
gated by iBulik et al.l 01 998) using the potentials calculated by 
iKuiiken & Gilmord (Il989,) . 



2.3. Lens search algorithm 

A naive search for lenses i.e. testing for events when a star and 
a compact object are aligned with Earth on its galactic orbit 
would require nearly infinite computational resources. Therefore 
we use the following lens search algorithm. We look for the lens 
events that take place anywhere on the galactic orbit of the Earth. 
This is much easier computationally as this approach make use 
of the cylindrical symmetry of the Milky Way. 

We assume that the Earth's orbit is circular with the radius 
of 8.5 kpc. Let us consider a sphere of radius R - ^.5 kpc (which 
is equal to Earth's galactic orbit radius) and centered in the cen- 
ter of the galaxy. Let us also consider a star-lens system located 
somewere in the galaxy. The position of the system is charac- 
terized by two vectors: P - the position vector of the star in the 
galaxy, and D - the star-lens vector (see Fig.lTJ. We search for a 
point Q in which the line passing through the position of the star 
and the lens pierces the sphere. In principle we look for a factor 
t which satisfies the condition 



\P + tD\=R. 
Eq. (IDi is a quadratic equation in t 
(P + tD) ■ (P + tD) = R^ , 



D^t^ +2D Pt + (P^ -R^)^0. 



(4) 

(5) 
(6) 
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Now, for each pair of a star and a compact object, we cal- 
culate A taking into account their positions at the beginning and 
at the end of a calculation step. If A > we calculate A' from 
the position of the pair at the end of the step, which in our cal- 
culations was one month. If it is also positive, then we calculate 
both fi_2 and f'j respectively. We are only interested in a con- 
figuration for which fi 2 > 1 and f j , > 1. If f < the compact 
object is located behind the star, and for < f < 1 the compact 
object and the star reside on the opposite sides of the sphere. In 
both such cases the lensing event cannot be observed on Earth. 
For both f and f' we then calculate positions of the points where 
star-lens line crosses the sphere at the beginning and at the end 
of the examined evolution period. The lensing event happens if 
these points are located in different hemispheres {z'z < 0). 

This algorithm is not sensitive to some special configurations 
when during the evolution star or compact object are located 
outside the sphere at very high azimuthal angles (as seen from 
Earth). We checked for these special occasions with two addi- 
tional algorithms. We conclude that such events are extremely 
rare (there was not even one such case) and does not influence 
our results. 



2.4. Lensing rate 

In order to estimate the true lensing rate one has to scale the 
number of the lensing events found in the simulation to match 
the physical conditions in our Galaxy. There are two scaling fac- 
tors involved. One results from the number of simulated stars 
and compact objects and the second from the fact that we look 
for the lensing events that happen on the whole Earth's orbit. 
Let us denote the Earth orbital period in the Milky way as 
P - 250Myrs.Let the true number of star in the Milky Way 
be A^*j X 10", and fie„s - represent the fraction 

of the stars in the galaxy that can be detected by lens search 
experiments, and the number of compact objects in the Galaxy 



A^^'^ X 10**. Our simulations include lO*" stars that can be lensed 
in the Galactic potential. 

Thus the expected rate of lenses due to the compact objects 

is 



R^4x 10^/;, 



leiis,- 



-1 V elm 



(8) 



where Nsim is the number of compact objects (potential lenses) 
in the simulation, and Niens is the number of lenses we find in 
one yean 

3. Results 

Using Startrack population synthesis code we evolved 10^ bi- 
nary systems, obtaining data about compact objects (neutron 
stars and black holes). We put those objects (only solitary, from 
disrupted systems) into the model Galaxy. We also randomly in- 
troduced 10^ stars in the Milky Way. All the stars were then 
evolved for IQGy in Galactic potential taking into account the 
initial velocities that the compact objects received at birth due to 
asymm etric kicks and disruptions of the binary systems Blaau^ 
( Il96ll) . We then searched for lensing events that take place any- 
where on the Earth orbit in the galaxy in one year of simulation. 
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Fig. 2. Skymap of the sumlated lensing events for the standard 
model of compact object formation. The skymaps in the two al- 
ternative models look similar. 

Table 2. The results of the lens search simulations and the ex- 
pected galactic lensing rates. 

Model 





^ ' sun 


^ lens 


Rate[yr'] 


A 


76399 


826 


4.3 


W0.5 


73575 


757 


4.1 


KO 


26619 


584 


8.7 


Kl 


73907 


687 


3.7 


C 


81327 


887 


4.3 


S 


60340 


1368 


9.0 


M 


8170 


194 


9.5 



The resulting rates, calculated using equation[8] are presented in 
Table 1 for each model assuming that A^*j = l,andA^^'' = l.The 
typical values obtained are a few tens per year, however these ex- 
pected rates have to be considered as rough estimates given the 
number of assumptions used in equation [8] Note that the Table 
gives the values of the rates assuming that all compact objects 
are formed in the particular scenario. 

In each simulation we have noted the positions of the lens- 
ing events on the sky. We present the skymaps of these posi- 
tion in Figure|2l The lenses are strongly concentrated around the 
galactic center. In Figure [3] we present the cumulative distribu- 
tion of the fraction of events as a function of the distance from 
the Galactic Center. Typically 70% of the lensing events take 
place within x 5° from the Galactic Center, while a fraction 0.9 
of the events happens within 20°. 

We also note the mass of the gravitational lens for each lens- 
ing event. In Figures|4]and|5]we present the intrinsic distributions 
of masses of compact objects for each model and the mass dis- 
tribution of lenses. In the case of models S, M, KO, and Kl the 
distributions are either quite similar. They do differ in the case 
of remaining models. But if they differ then they do in the same 
manner:the distribution of observed masses shows a deficiency 
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Fig. 3. Cumulative distributions of the angular distance lensing 
events from the Galactic Center for the models considered in the 
paper. 
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Fig. 4. The cumulative distributions of masses of compact ob- 
jects in the standard model A. The thick line denotes the intrinsic 
distribution of compact object masses while the thin line corre- 
sponds to the distribution of measured masses of the lenses, each 
plot is labeled by the model as in Table 1 . 



of low mass compact objects (neutron stars) and an increased 
fraction of high mass objects (black holes). The distributions are 
similar for the models in which the velocities of compact objects 
weakly depend on their masses: KO, and Kl. 

In the cases when the distributions differ the suppresion of 
low mass compact objects is not more than about 20% in relation 
to the high mass compact objects. 




Fig. 5. Same as Figure 4 but for the remaining models. 



4. Summary 

We have considered several scenarios leading to formation of 
galactic compact objects and calculated their initial velocity dis- 
tributions. We then simulate the spatial and velocity distributions 
of them in our Galaxy. Each model yields a different compact 
objects mass distribution and may have a different distribution 
of velocities. We simulate the microlensing events and note both 
mass and position for each one. 

We estimate the rates of lensing events due to compact ob- 
jects. These estimates show that they can be observed. The rates 
in Table 2 must be taken as very rough estimates. There are some 
poorly known factors that are included in these rates - see equa- 
tion [8] One should note that the fraction of objects detectable 
by the lens searches //f„j can be increased for the searches con- 
ducted in the infrared. However, we know that there already are 
some events which can be interpreted as microlensing events by 
black holes. 

In the simulations the robust result is the distribution of lens- 
ing events in sky. All the events are well concentrated around 
the Galactic center Thus the lens searches should concentrate 
on this region. This is probably due to the fact that the dominant 
factor determining the position on the sky is the density of the 
lensed stars and not lenses themselves. 

Finally, we find that the observed mass distributions differ by 
less than 20% from the intrinsic one. The observations of several 
gravitational lenses from black holes or neutron stars should lead 
to relatively good estimate of the intrinsic compact object mass 
distribution. The observation of the mass distribution of compact 
objects is extremely important since this is a unique way to probe 
the final stages of the massive stars evolution. Observations of 



black holes in massive accreting binaries may suffer from nu- 
merous selection effects and probably do not reflect the underly- 
ing intrinsic distribution. 

The lensing events due to massive lenses will have long du- 
ration and therefore require long observational campaigns. We 
emphasize the importance of such searches, as their potential 
results will have a very significant relevance for the theory of 
compact object formation and models of massive star evolution. 
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